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Abstract 

For parabolic spatially discrete equations, we consider Green’s functions, also known 
as heat kernels on lattices. We obtain their asymptotic expansions with respect to pow¬ 
ers of time variable t up to an arbitrary order and estimate the remainders uniformly 
on the whole lattice. The spatially discrete (difference) operators under consideration 
are finite-difference approximations of continuous strongly elliptic differential opera¬ 
tors (with constant coefficients) of arbitrary even order in M. d with arbitrary d € N. 

This genericity, besides numerical and deterministic lattice-dynamics applications, al¬ 
lows one to obtain higher-order asymptotics of transition probability functions for 
continuous-time random walks on Z d and other lattices. 
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1 Introduction 

The paper deals with Green’s functions, or heat kernels, of general parabolic equations with 
constant coefficients that are continuous in time and discrete in space. Applications that we 
have in mind include spatial discretizations of continuous models, lattice dynamical systems, 
and continuous-time random walks. 

We consider parabolic problems on the grid space , or lattice, 

Mg := {x G R d : Xk = Sk£, G Z, k — 1,..., d}, d e N, e > 0, 

of the form 

u £ (x, t) + A e u £ (x, t) = 0, x G Mg, t > 0, 
u £ (x, 0) = 5 £ (x), x G Mg. 

Here ' = d/dt, 8 £ (x) is the grid delta-function given by 

<5 £ (0) = £ ~ d , 8 £ {x) =0 Vx G R d e \ {0}. 

and A e is an elliptic difference operator (with constant coefficients) which is assumed to be 
an Mth order approximation of a strongly elliptic (continuous) differential operator A(T>) 
of even order £ (see Section [2] for rigorous definitions). An order M 6 N is defined by the 
estimate 

sup \A e u(y) — A(V)u(y)\ < C{u)e M \/e > 0, (1.3) 

y£R d 

which must hold for any smooth (continuous and bounded with all its derivatives) function 
u(y) and an appropriate constant C(u) > 0 not depending on e > 0. 

We call a solution of problem (11.111 the first discrete Green function. Given arbitrary 
initial data instead of 8 s fx), one can represent a solution as a (discrete) convolution of 
these initial data with the first Green function. Alternatively, the first Green function 


( 1 . 1 ) 

( 1 . 2 ) 
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determines a semigroup generated by — A e . Various properties of A e (in terms of its symbol) 
in functional spaces of grid functions, corresponding semigroup properties, and relations to 
spatially continuous problems can be found in mn, see also references therein. 

By using the discrete Fourier transform, one can give an integral representation of the 
discrete Green function (see H2.2211 ). However, it is not possible to express it via elementary 
functions, hence its asymptotic expansions play an important role. Asymptotics of the so- 
called lattice Green’s functions in the stationary case were studied beginning from 1950s, 
see [5], the subsequent papers and the monograph ns Chapter 8] . For 

parabolic operators, there is vast literature in the spatially continuous case. For example, 
large-time behavior of Green’s functions was treated in [18] (for small perturbations of the 
heat operator) and in [19jl22j (for spatially periodic coefficients). A survey on the large time 
behavior of heat kernels for second-order parabolic operators on Riemannian manifolds can 
be found in [2T] . In the spatially discrete case, the research directions include continuous¬ 
time random walks on general graphs (see, e.g., [I3li20] and references therein) and on lattices 
in a random environment (see, e.g., f4j). In both cases, Gaussian bounds for the heat kernel 
is an important question. However, higher-order asymptotics of Green’s functions is not 
available in general. We mention men] , where an asymptotic expansion of Green’s function 
for specific parabolic equations on one-dimensional lattices was obtained in terms of the 
Bessel functions. 

In our paper, using the integral Fourier representation of the first Green function u £ (x,f), 
we obtain a higher-order asymptotic formula of the form 

<9 J u e (x, t) 1 / x \ 

frl = JdjM Hj {fUl) 

K k (1.4) 

+{k+d)/i+J H Jk (fTT?) r u(^’ H,X,t), x £ ®e> t — hF • 

k=M 


Here £ > 0, to > 0, K > M, J > 0, Hj(y ) and Hj k (y), y G M d , are explicitly given smooth 
functions, and the remainder satisfies 


r £ u (J,K;x,t)\ < 


£ K+ 1 R u (J,A,t 0 ) 

t (K+d+i)/e+j ’ 


x G Mf, t > t 0 e e , 


(1.5) 


with R U (J, K, t 0 ) > 0 not depending on £ > 0, x G and t > to£ e . We emphasize that the 


estimate in (II. 5ft is uniform with respect to x G 

Moreover, we show that, for J = 0, the leading order term 


1 

t d/i 


H 0 


in (II.4p coincides 


. (£) 

with the continuous parabolic Green function (see (I2.6p h Hence, (1 1.4 p implies that the 
first discrete Green function approximates the continuous one, uniformly for x G Mf, with 

an error O (^ ^ M+d y^ j as t —> oo and with an error 0 (£ M ) as £ —> 0, where M is the 

order of approximation in (11.31) (see Corollaries 13.11 and 13.21) . In the case where — A e is an 


approximation of the Laplacian, the first term in (11.41) with J = 0 equals 


1 


{2\/nty 


_ e —M 2 /4i_ 
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In particular, it is radially symmetric. The other terms, being non-symmetric, take into 
account the radial non-symmetry of the lattice Rf. 

Asymptotic formula (jl.4j) is proved in Section [HI (Theorem 13. l|i . after rigorous definitions 
of elliptic differential and difference operators in Section [2j 

We note that the general form of the difference operator A E (rather than the Laplacian 
only) allows one to treat continuous-time random walks on d- dimensional lattices that are not 
necessarily cubic, but still represented by a discrete additive subgroup of R d . This is possible 
whenever one can linearly transform the vertices of such a lattice to R d (see [13, Section 1]) 
and obtain an elliptic generator A £ on Rf. Then, linearly transforming (II.4ft back to the 
original non-cubic lattice, one can deduce asymptotics of the transition probability function 
of the original random walk. Example 12.11 2 illustrates an elliptic difference operator that is 
a generator obtained after transforming the random walk on the triangular lattice to 71?. 

Along with the first discrete Green function, we study the second discrete Green function 
in Sections [4] and [5j It is defined as a solution of the problem 

f v £ (x, t) + A e v(x, t) = S £ (x), x e Mg, t > °> gs 

1 v £ (a:, 0) = 0, x E R d . 


Given an arbitrary time-independent right-hand side instead of S £ (x), one can represent a 
solution as a (discrete) convolution of this right-hand side with the second Green function. 
An interesting application of the second Green function v £ (x,t) occurs in discrete reaction- 
diffusion equations with hysteretic nonlinearity in the right-hand side. Higher-order asymp¬ 
totic formulas for v £ (x, t ) play a central role in understanding pattern formation mechanisms 
there (see 0 )- 

Asymptotics of the second Green function depends on the sign of £—d, where £ is the order 
of the continuous differential operator AifD) and d is the spatial dimension. If £ — d < — 1, 
then we prove that the second discrete Green function satisfies 

V'M) = Fo (^) + phsj(£) 

K k (1.7) 

+{k+d)/i-l ( /1 /£ ) r v(-^) x A)i x e t > to£ , 

k=M 


where the remainder is estimated as 

K(K] x,t) | < £ t(K + d l[ye-i ^ » x eR d £ , t>t 0 e l , ( 1 . 8 ) 

with R v (K,to) A 0 not depending on £ > 0, x G M^, and t > t 0 £ e . Again, the estimate 
in (II.5p is uniform with respect to x G Mf. In Section |4] (Lemma 14.21) , we prove (JTTTJ) with 
the functions F k (y) that are given explicitly as solutions of certain first-order PDEs with 
H 0 (y ) and H ok (y) in the right-hand side, but still unknown function In Section 15711 

(Theorem 15.ip . we find the function ?l(y) by simultaneously passing to the limit, as t —> oo, 
in ([L7l) and in the explicit integral representation of v £ (x,t). 
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If £ — d > 0, the situation is more delicate. In Section 0 we prove an analogue of (jLTJ) 
in the form 


e-d 


v'OM) = (0) + ]T eV-™/' Ft (0) + £ '- d n 0 

(JL\ 


k=M 
k 


K £ k (1 ' 9) 

+ Y M+im -1 Ft iw) + ^ K ' x ’ *)’ 1 e R ? \ {0}. t > 

k=max(M/-d+l) 


(if M > l—d+1, the first sum is absent). As before, the remainder satisfies (j 1.8 [) . the functions 
F k (y) are given explicitly as solutions of certain first-order PDEs with H 0 (y ) and H ok (y ) in 
the right-hand side, and Yl(y) is still unknown. However, unlike in the case £ — d < —1, 
formula (1 1.9 j) may contain the function F k -d{y), which in general appears to be undefined at 
y = 0 and to have logarithmic growth as y —> 0. Furthermore, we cannot pass to the limit as 
t — > oo immediately because both the terms with k — 0, M,... — d in (11,9p and the explicit 

integral representation of v e (x, t) tend to infinity. Hence, in Section [572] ('Theorem [572]) , we use 
the explicit integral representation of v £ (x,t) to deduce another asymptotic representation. 
It contains a linear combination of nonnegative powers 0 — 0 — d) and of 

Inf with ^-dependent coefficients and holds for each fixed (but not uniformly in M^). 

Comparing it with (j 1.9 [) allows one to determine Q(y) and obtain uniform asymptotics in 
In Section [6l we apply our general results to the (27V)th order approximation of the 
one-dimensional Laplacian: 

N 

A £ = —A £ , A £ u(x) := e~ 2 a !/ (u(a; — ev) — 2u(x) + u(x + ev)) (1.10) 

U=1 


with appropriately chosen coefficients a u G R. (In this case, £ = 2 and d — 1, hence we are 
in the situation £ — d > 0.) First, we claim that A e is elliptic for any N G N in the sense of 
Condition 12.31 (the proof is given in Appendix O. Then we explicitly find all the functions in 
expansions (jl.4[) and (j 1.9 [) for the first and second Green functions and in the corresponding 
expansions for the (spatial) gradients of the first and second Green functions. 

Interestingly, if N = 1, i.e., A e u(a:) := £~ 2 (u {x — e) — 2u(x) + u(x + e)), it appears 
that VL(x/e) in (11.911 vanishes for all x G Rf, which is proved in Appendix [B] In general, 
VL{x/e) ^ 0 for N > 2. 

To conclude, we collect some general notation that we use throughout the paper. We use 
variables x G R £ or x G 7L d (for e — 1) and y G R d for the grid and continuous coordinates, 
respectively, and rj G R ne - 1 , 0 G R n , and (Gl d for the Fourier coordinates, where 

R x :={£eR d : |e fc | < A, fc = l,...,d}, A > 0. (1.11) 

By B\, we denote the ball of radius A > 0 in R d . We denote by (r, (pi,... (pd-i) the spherical 
coordinates of y G R d or 6 G R n and introduce the function 


r n (<p) := distance from the origin to the boundary of R^ in the direction (p. (1.12) 
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We use the upper-case bold letters for difference operators: A e , A_- and lower-case bold 
letters for grid functions (i.e., functions defined on M^): u £ (a;,f), v £ (x,t), etc., indicating by 
superscript £ that x G Rf. We may omit the superscript if e = 1 and x G 7L d . We denote 
polynomials by calligraphic letters: A(£), V(£), etc. and other (non-polynomial) functions 
by usual letters: A{6), H(y), fl(y), u(y,t), etc. 


2 Continuous and discrete parabolic Green’s functions 

In this section, we introduce parabolic equations with general elliptic operators of even 
order t in R d , d G N, and the corresponding finite-difference approximation and define their 
Green functions. 


2.1 Continuous Green function 

Consider a differential operator with constant coefficients in R d , d G N: 

A(V) = Y, b a V a , (2.1) 

\a\=l 

where £ G N is even, a = (aq,. .., otd) is a multi-index, \a\ = aq + • • • + a^, b a G C, 
V a = Vf 1 ... Vf d , and Vj = -id/dyj, j = 1,..., d. 

Definition 2.1. The polynomial 

■4({) := J2 e 6 < (2.2) 

\a\=i 

where .(obtained from A{T>) by formally replacing Vj by fj), is called the 

symbol of the differential operator A{T>). 

Obviously, there is C > 0 such that 

|A(0|<C|^ VfeR d . (2.3) 

We assume throughout that the operator A(V) is strongly elliptic, i.e., the following 
condition holds. 


Condition 2.1. There is c > 0 such that 

Fte-4(0 > e|0' V£eR d . 


(2.4) 


We briefly recall the notion of a continuous Green function for parabolic equations. We 
define it as a solution of the problem 


du(y,t) 
dt 

u(y, 0) = S(y), 


+ A(V)u(y, t) = 0, y G K d , t > 0, 

yeR d , 


(2.5) 
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with the 5-function in the initial condition. To give an explicit formula for the Green function, 
we formally use the Fourier transform and deduce from (12.51) 


u(y,t) = —^ [ e tA{ ®e w<i d£, y eR d , t > 0, 

G 7r J J R d 

where the integral converges for each fixed t > 0 due to (12.41) . Making the change of variables 
£ H y and using the homogeneity of M(£), we obtain 

u (v, t) = -^i H j t ) > yeR d ,t> 0, (2.6) 


where 

H (V)-=J .dw / (2,7) 

(2vr)“ J R d 

Due to (I2.4|) . H(y) and all its derivatives decay at infinity faster than any negative power 

of \y\. 


Definition 2.2. We call the function given by (12.6[) the continuous Green function. 


Remark 2.1. If A(T>) = —A, where A is the Laplace operator, then M(£) = |£| 2 , H(y) = 

———e — l^l 2 / 4 , and (12.6p assumes the well-known form u(y,t ) =--==—e - ^ 2 / 4 *. 

2 d 7 T d /- (2 \m) d 


2.2 Discrete Green function 

For e > 0, we define the grid space 

:= {x G R d : x k = s k £, s^eZ, k — 1,..., d}. 

A function defined on R d is called a grid function. We say that a grid function (continuous 
function) is rapidly decreasing if it decays at infinity faster than any negative power of |x|, 
igRf (of \y\, y G R d ). For a continuous function u(y), y G R d , we use the notation 

8 e ,k±u(y) = u(y ± e k e) - u(y), (2.8) 


where e k is the unit vector of the axis y k . The same notation is used for grid functions u £ (x), 
x G Rg. 

We fix a natural i > l and consider the following difference operator with constant 
coefficients: 


A e := 




^ 2 d-l cv 2 d 
°£,d~ °£,d +’ 


(2.9) 


i<\v\<i 


where v = (ui, ..., u 2 d) is a multi-index, \v\ = + ■ ■ ■ + u 2 d, and a u G C. 

We say that a function u(y), y G R d , is smooth if it is continuous and bounded together 
with all its derivatives. 
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Definition 2.3. For M £ N, the difference operator A £ is said to be an Mth order approxi¬ 
mation of the differential operator A(T>) if, for any smooth function u(y), there is a constant 
C{u) >0 such that 


sup |A e u(y) — A(V)u(y)\ < C{u)e M We > 0. (2.10) 

y£ K d 

We assume throughout that the following holds. 

Condition 2.2. There is M £ N such that A £ is an Mth order approximation of the 
differential operator AifD). 

Next, we introduce a symbol of A £ and relate it to the symbol of AifD). To motivate 
the notion of a symbol, we define the Fourier transform of a rapidly decreasing grid function 
u e (x) by the formula 


(F £ u £ )(^):=(27r)- d ^4-*, (2.11) 

where R ^ e -1 is given by fll.lll) . The Fourier transform F £ establishes an isomorphism between 
rapidly decreasing grid functions and smooth functions that are 27re _1 -periodic with respect 
to each of the variables rji, ..., yd- The inverse transform is given by 

(F" 1 ^)^) = f e lxr, v{rj) dy. (2.12) 

If u £ (a;) is a rapidly decreasing grid function, then A e u e (a;) is also rapidly decreasing and 

(F £ A £ u £ )(t?) = £~ e A(ye) ■ (F e u £ )(y), (2.13) 

where the function £~ l A{ye) is obtained by replacing the operator S Et k± in the right-hand 
side of equality (12.9h by the expression e ±ir,kE — 1. 

Definition 2.4. The function 

A(0):= a » {e~ i61 ~ I)" 1 {e iei - I)" 2 - - - (e~ i9d - l)" 211 - 1 (e i6d - l) U2d , 6 £ R n , (2.14) 

e<\v\<£ 

is called the symbol of the difference operator A £ . 

The following lemma establishes a connection between the symbols of AfD) and A e . 
Lemma 2.1. If Condition al . 21 holds, then, for any K > M, 

K 

A{e) = Y,A t+ m + o(\e\ K+l+1 ), 

k =0 
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(2.15) 








where 



r m 

CD 



Ak+i{0) < 

0 

ifk = 1, ■ ■ 

., M — 1, 

(2.16) 

homogeneous polynomial 


of degree k + £ 

if k — M,. 




and O(-) is understood with respect to \6\ —> 0. 

If the symbol of A e satisfies (12.15[) and (12.16[) with some K > M, then Condition 12.21 
holds. 


Proof. We observe that, for any smooth u(y) and any J > 1, 


■w<(») = E 

e*"* - 1 = E + O(|0| J+1 ), 


i=i 


r- 


where sup \Rk,j+i(y)\ < C'i(w) with some Ci(w) > 0. Therefore, taking into account defini- 

j/gR d 

tions (12.91) and (12.141) . we have for any K > M 


K 


A eu(y) = ^2£ k A k+ e(V)u(y) + e K+1 R K+e+1 (y), 

k =0 
K 

A(e) = J 2 ^i( 0 ) + o(\e\ K+w ), 


k =0 


where Ak+e(-) are the same homogeneous polynomials of degree k + £ in both formulas and 
sup \RK+e+i(y)\ < C 2 (u) with some C 2 (u) > 0. Comparing the last two equalities and 

ygR d 

taking into account (j2. 1Q|) . we conclude the proof. □ 

Along with Condition 12.21 we assume throughout that the following ellipticity condition 
holds. 

Condition 2.3. Re A{6) > 0 for all 6 e R n \ {0}, where R n is defined in (11.111) . 

Remark 2.2. Condition 12.31 does not automatically follow from the fact that A e is an 
approximation of a strongly elliptic operator A(T>) (only the inequality Re A(6) > 0 for 
small 1 9 1 ^ 0 follows). However, Lemma 16.21 below implies that Condition 12.31 holds, for 
example, for any order approximation of the second derivative in R (see (j6.8j) and (16.9jl ) and, 
hence, for the corresponding approximation of the Laplace operator in R d , d G N. 
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Example 2.1. 1. Consider the difference operator 

d 

A e := £ y~^ 5 e j+6 e j-. 

3 = 1 

For e = 1, the operator —(2d) _ 1 A! is a generator of the continuous-time random walk on 7L d 
(see [121 Section 1]). The operator A £ is the 2nd order approximation of the strongly elliptic 
operator A(T>) = —A, where A is the Laplacian. The symbol of A e is given by 


A (9) := 2^(1 - cos Oj). 

3=1 

In Section [ 6 l we will also consider (2N)th order approximations of the Laplacian for d — 1. 
2. Let d — 2. Consider the operator 


A« “ ^5(24,1+4,1- + 24,2+4,2- - 4,1+4,2- - 4,2+4,!-). 


(2.17) 


The operator A, is the 2nd order approximation of the strongly elliptic operator A(T>)u : = 
4 

~-^( u yiyi + u yiy 2 * Uyiyi)- Tlie symbol of A e is given by 

A{9) — 4 — - (cos 6\ + cos 6*2 + cos($i — 0 2 )) 
o 

and the symbol of A(T>) is given by 


A(d) = -(9f + 9l + (9 1 ~9 2 y). 

On the other hand, — A e corresponds to the discretization on the triangular lattice of the 
continuous Laplacian. Indeed, the matrix 


M : = 


1 -l/\/3\ 

0 2/^3 ) ' 


takes the triangular lattice to 7? and, at the same time, A((M T ) : £) = The latter 

is the symbol of the differential operator —A. 

Lemma 12.11 together with Conditions 12.1112.21 and 12.31 implies the following. 

Lemma 2.2. We have 

|A(0)| < C\9f V9 e R w , (2.18) 

ReA(9)>c\9\ e V9 e R n , (2.19) 

where C and c are positive constants that do not depend on 9 <E R^. 
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Without loss of generality, we assume that the constants C and c in (12.18jl and (12.1 9ft are 
the same as in (12.3p and (12.41) . respectively. 

Now we introduce discrete Green functions. Let 8 £ fx) be the grid delta-function , see (jl,2j) . 

Definition 2.5. We call the function u £ fx,t), x G Mg, t > 0, th e first discrete Green function 
if u £ (-, t) is a rapidly decreasing grid function for all t > 0, u £ (x, •) G C' 1 [0, oo) for all x G Mg, 
and 


u £ (x, t ) + AgU £ (x, t) = 0, x G M(, t > 0, 


u e (x, 0) = S £ (x), 


x G 


( 2 . 20 ) 


We call the function v £ (x, f), x G Mg, t > 0, the second discrete Green function if v £ (-,t) 
is a rapidly decreasing grid function for all t > 0, v £ (x, •) G C' 1 [0, oo) for all x G Mf, and 


v £ (x, t ) + A e v £ (x, t ) = (5 £ (x), x G M(?, t > 0, 
v £ (x, 0) = 0, x G Mg. 


( 2 . 21 ) 


Here' = 9/cif. 


In Section [3l we will establish a connection between the discrete and continuous Green 
functions in terms of asymptotic formulas. 

Using the discrete Fourier transform (12. lip , its inverse (12. 12p . and relation (12. 13p . we 
obtain the explicit representations 


u £ (x, t ) = v £ (x, t ) = 

1 f 


IR _i 

7T £ - 1 - 


v £ (x,t) = 


(2tt) 


(2n) d 
r x — e - te ~ l M r i £ ) 
Ir _i £~ l A{r]£) 


e -ts e A( V ) e ix V dr} , 


e ixri d V . 


Changing the variables in the integrals in (12.221) and (I2.23p . we obtain for J 


( 2 . 22 ) 

(2.23) 

0,1,2,... 


<9 J u e (x, t) _ _ J£ _ d d J u 1 (x / ,r) 

dt 1 =£ d7 J 


1 

x'=x/e,T=t/e l 


V £ (x, t ) = £ f ^V 1 




(2.24) 


3 Asymptotics of the first Green function u £ (x,t) 


3.1 Formulation of the result 


In the theorems below, we will not explicitly indicate the dependence of the remainders in 
asymptotic formulas on the number M fixed in Condition 12.21 

In this section, we obtain the following asymptotics for u £ (x,f). 

Theorem 3.1. For any £ > 0, to > 0, integer K > M, integer J > 0, and all x G Mg and 
t > t 0 £y we have 


d J u £ (x, t ) 
dt J 


E 


t (k+d)/e+J 


H Jk 



+ r £ {J,K\x,t), 


(3.1) 
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where 


H Jk {y) < 

H(y) is given by (12.711 . 


((-A(V)) J H(y) ifk = 0, 

0 if k = 1,, M — 1, 

finite linear combinations 

of derivatives of H(y ) if k = M ,..., AT, 


(3.2) 


\r e u (J,K;x,t)\ < 


e K+1 R u {J,K,t 0 ) 


(3.3) 


£(tf+d+l)/£fJ 

and R U (J , A", t 0 ) > 0 does not depend on £ > 0, x G , and t > t 0 £ e . 

Remark 3.1. It follows from the proof below that Hjk(y ) = TZj^+e('R)H(y), where 
are polynomials of degree >k + l defined by (13.15(1 . 

Note that the main term in the asymptotic representation of u e (x,t) coincides with the 
continuous Green function (12.61) . Theorem 13.11 implies that the higher the approximation 
order M (see Definition 12.3j) is, the faster the first discrete Green function converges to the 
continuous Green function, as t — > oo or as £ —» 0. More precisely, the following corollaries 
hold. 

Corollary 3.1. Under the assumptions of Theorem 13.11 for any £q > 0, 

±h(J£\ 


sup sup 

£G(0,£ 0 ] zGKf 


u ( x fi) yd/1 H \u/i J 


= 0 ( t (M+d)/i ) ast^oo. 


Corollary 3.2. Under the assumptions of Theorem 13.11 for any t\ > 0, 


sup sup 

t>ti xGRd 




= 0(e M ) as e —> 0. 


3.2 Proof of Theorem CUE 

Setp 1. First, we assume that s —— 1 and denote ui— Note that |n (x,t)\ 

is bounded for x £ 7L d and t > 0 due to (12.221) and (12.191) . Therefore, if we show that 
Theorem 16.11 holds with £ — 1 and some to > 0, it will follow that the theorem holds with 
£ — 1 and all to > 0. Let us fix t 0 = 1 and consider t > 1. 

Differentiating (12.221) with £ = 1 and making the change of variables £ = t l ^rj 1 we write 

d ~^r 1 = WfW / { ~ 1)J W rl/ '«)' - eixr '" ( =« x ’ ( >+^ *>, 

(3.4) 

where 

/l(l ' i);= (2 WW / (-1 Y di, (3.5) 
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Uz.t) ■■= (2w) ] tJ/t J (~1) J (Ait-^oY (3.6) 

R *tV*-\ R *tV L 

L:=(£ + 2)1. (3.7) 

Step 2. Let us estimate h(x,t ) in (13.511 . First, we consider the term e tA! ' t By- 

Lemma 12.11 

K 

= A(® + y t-‘/'A + ,K) + i(i, ®, (3.8) 

k=M 

where 

\AM\ < Cl r^ +1 )/ < |e| K+1+< , £ e t > l. (3.9) 

Here and below c±, C 2 ,... > 0 do not depending on £ and t. 

Further, (13.7ft implies for £ G R nt i/z, and t > 1 that 

< c ^ t -k/z t (k+i)/L = c ^(£-(e + i)k)/L k = 1,2,.... 


Therefore, 


K 


yt- h ' ,c A M (o+ 


k=M 


<c 3 , ( e < > l. 


(3.10) 


Using the Taylor expansion of the exponent function, we obtain 


K 


K 


K 


k=M 


exp ( - y - Ait.® ]=l+^ s n t- k/, A t+ ,(i) + i(0 ) + ilK.i) 

m= 1 ' \k=M 

K 

= 1 +Y t t~ k/e ' p k+i(0 + P(U) 


k=M 


(3.11) 

Here Vk+e(Q are polynomials of degree > k + £ and the remainders are estimated due to 
relations (13.9ft and (13.101) as follows: 


\MZ,t)\ <e C3 


K 


yt- k/c A M d)+A(i,t) 


k=M 


K +1 


\P(Z,t)\ < cse^'H-V+W, £ G R nt 1/,, t > 1, 

where c is the constant from (12.41) . 

Now we estimate the term (H(t _1 ^£)) J . Using (13.81) , we have 


< c,e c ^/ 2 t- M ^ K+1 ^, £ G R nt i/L, t> 1, 

(3.12) 


K 


{Air 1 ! 1 ®) = r J MK)) J + £ r h/l Vj, h+ e(® + D(t,® , 


(3.13) 


k=M 
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where Vj^+eiO are polynomials of degree > k + £ and 

\D(U )I < cee^t-V+W £ e R^ /h t > 1. 

Combining (liUTD - fliUlD with (MJ) yields 

(-1 y (Air^oY e- tA{t ~ y ^ 

f K 

= e- A[i) ( t- J (~A(0) J + Y, + R(t,a I, 


(3.14) 


(3.15) 


k=M 


where 7 Zj^+e(0 are polynomials of degree > k + £ and 

m,t)\ < c 7 e c ^/ 2 t-^ K+i ^- j } e e r^i/l, t > i. 

Substituting (13.151) into (13.5|) . we have 

1 r (—^K)) J e-- 4 < { ) ■e ixr ' l ‘ i di 


(3.16) 


I\ (x,t) = 


(2Tc) d t d / £+J 


R 


t l/L 


K 

E 

k=M 


(27[) d ^~\~J 


Kj, t+ ,K)e^ Ai) -e ixrl '‘ ( d4 ( 3 . 17 ) 


R 


(2'K) d t d R 


>*!/•£- 

e~ A ^R(C,t) -e ixt ~ 1/l t dC 


R ntl It 


Note that, due to (12.4[) . the polynomials £?(£) = (—4(£)) J and £>(£) = IZj^+iiC) satisfy 


■ e ixt ^ d£ 


t l/L 


< 


\m\ ■e-W t d£ = 0 




t (K+i)/e I ’ 


(3.18) 


where O(-) is taken as t —> oo (uniformly with respect to x ). Similarly, due to (13.16j) . R(£,,t) 
satisfies 


e~ A ^R^,t) -e^'^di 


R 


t l/L 


< 


< 


C 7 


t (K+l)/l+J 


3 -Re-4(f) . _c|€|V2 


e c|41 /z d£ 


\tV L 


(3.19) 


Si _ / ,-d€l72 


t {K+l)/i+J 


d£. 
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Combining fj3.17jl - fl3.19p . we obtain 


Il(yX '^ (2n) d t d /^ J 


J 





(3.20) 


Relations f!2.7p and f!3.20p and the fact that 1Zj^+z{C) are polynomials imply 



K 


1 


t {K+d+i)/e+j 


(3.21) 


k=M 


where H J0 (y) = (-A(V)) J H(y) and H Jk (y) = 1Z J}k+i ( V)H(y ). 

Step 3. Let us estimate / 2 (x,t) in ( I3.6j) . Using (12.18)) and (I2.19p . we have for £ G 


\ ^nt 1 / L 


(Air^oY e- tA[t ~ 1/e ° < c 8 e~ c < c 8 e-^* IL l 2 


Therefore, 


I 2 (x,t)\ < c 9 e- e * t *' L / 2 , t> 1 . 


(3.22) 


Relations fl3.4p - fl3.6p . f!3.2ip and f!3.22p yield the assertion of Theorem 13.11 with e — 1. 
Using the first relation in (I2.24p . we obtain the assertion with any e > 0. 

4 Preliminary asymptotics of the second Green func¬ 
tion v £ (x, t) 

Throughout this section, we assume that £ = 1 and denote v(x,f) := v 1 (x,t). First, we 
obtain an asymptotic formula for v(x,t) up to some unknown function fl(x). In the next 
section, we will provide an explicit formula for it. 

We introduce the spherical coordinates (r, p) with r > 0 and p £ 4>, where 


<f> := {p e R d 1 : pi, ■ ■ ■, <Pd- 2 e [o, 7 r], <p d - 1 e [0, 27 t)}. 
They are related to the Cartesian coordinates y e via 


(4.1) 


yi — r cos p \, i )2 — r sm ip i cos p 2 , 1/3 = r sm p\ sm ip 2 cos p 3 ,..., 

Vd-I = T' sin p 1 ■ ■ • sin p d _ 2 cos p d _ 1 , y d = r sin </?i • • • sin y 9 d _ 2 sin 


(4.2) 


and the Jacobian is given by 


dy = r d 1 J(p)drdp , J(</?) := sin d 2 </?isin d 3 y > 2 • • -sin^^. 


(4.3) 
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We will denote functions written in Cartesian and spherical coordinates by the same 
letter. For example H(r,<p) stands for H(y), etc. 

We set 

H k (y) := H ok (y ), y £ K d , 

where Ho k (y) are the functions from Theorem 13.II with J — 0. For y £ R d \{0}, we introduce 
the functions 

Ftiv) ■■= tr‘- d - k [°° dp ifk = 0 ,...,e-d, (4.4) 

Jr P 

U(s >)■■=—I'this,?) dp i!k>t-d+l. (4.5) 

r Jo 

Definition (14.41) will be used only in the case £ — d > 0 and definition (14.51) only for k > 0. 
Remark 4.1. Due to (13. 2p . F k (y) = 0 for k — 1,..., M — 1. 

Lemma 4.1. Let a function F k (y) be defined by (14.41) or (14.51) . Then 

1. it satisfies the partial differential equation 

d f)F 

(£-d- k)F k - Vi-Q 1 = iH k, !/eK'\{0}. (4.6) 

j =i yj 

2. it tends to 0 as \y\ —> oo, 

3. if k £ — d, it is bounded in W 1 and, for each fixed p, continuous with respect to r at 
the origin, 

4. if k > £ — d + 1, it is continuous with respect to y at the origin and 

U(0):=limU( !/ ) = - FT 4_^.(0). (4.7) 

Proof. Item Q] follows by observing that the equation in (14.61) is equivalent to the following: 

(£- d-k)F k -r-^- = £H k , r > 0. 

The latter can be solved as an ordinary differential equation for each fixed ip £ $. ft is 
easy to check that the unique solution that is bounded at infinity is F k (y) defined by (14.41) 
or (14.51) . respectively. 

Items [21 [3j and |4] follow from the fact that H k (y) are infinitely differentiable at the origin 
and are rapidly decreasing functions. □ 

Now we find a preliminary asymptotic formula for v(x,f) in terms of the above functions 
F k (y). ft contains an unknown grid function D(;r) and, if £ — d > 0, it is not valid at x — 0. 
The reason for the latter is essentially the fact that the functions F k (y ) are not continuous 
at y — 0 for k — 0 ,..., £ — d. 
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Lemma 4.2. There exists a grid function 12(x) such that, for any t 0 > 0, integer Ii > 
max(M, i — d), x G Z d \ {0} ; and t > t 0 , 


K 

v(x, t ) = 12(x) + 

k =0 


l 

t (k+d)/e -i 


F k 



(4.8) 


Here F k (y) are given by and (14.5p . 


r v (K;x,t)\ < 


Rv{K,t 0 ) 

t {K+d+l)/£-l ’ 


(4.9) 


and R V (K, t 0 ) > 0 does not depend on x G and t > to- 

Furthermore, if £ — d < —1, then the equality in (14. 8 p and the estimate in (14. 9 p also hold 
for x = 0 with F k ( 0) defined by (14.71) . 


Proof. We integrate (13. ip with J = 0 and e = 1 from t 0 to t and use the identity u(x, t ) = 


K f 

v(x, t) -v(x,t 0 ) = 

k =o ; 
to 


1 TT ( X 

T (k+d)/e k \ T i/e 


dr 


(4.10) 


+ / r^(0, K ; x, r) dr — / r^(0, K\ x, r) dr, x e Z d 


^0 


It follows from the chain rule and (14.6p that 

r (fc+d)/f ^ (tVl) = ( 4 ' X1 ) 

Combining (14. 10)1 and (14.111) and taking into account (13.31) and the inequality K > £ — d 
yields 


v(x,t) = Q(x) + t (k+d)/e-i Fk {£ji) + r vCftdM), xeZ d \{0}, (4.12) 

k =0 

with some grid function 12 (x) and the remainder r v (/l; x, t) satisfying (14 .9 [) . This proves (j4.8[) 
in the case x^0. If £ — d < — 1, the same formula for x = 0 follows directly from (14. 10)1 by 
substituting x = 0 therein and using (14.71) . 

To see that 12(x) and r v (K-,x,t) do not depend on t 0 , denote by 12* (x) and r *(K-,x,t) 
the functions corresponding to t 0 replaced by some t^ > t 0 . Then (14. 12jl implies that 


0 = 12(x) — 12* (x) + r V (K; x, t ) — r* (K ; x, t ), t > £g. 


Passing to the limit as t —» oo yields 12(x) = 12* (x). In turn, this and (I4.12p imply that 
r V (K; x, t ) = r* (K ; x, t ) for t > fg. □ 
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In the case £—d > 0, we will need modified versions of the functions in (14.41) . To introduce 
them, we denote by n^ the unit vector originated at 0 and pointing in the direction ip and by 
/bl (0, ip) the jth directional derivative of a function f(y) in the direction n v at the origin. 
Now, for each fixed ip € <3>, we expand Hk(p, <p) into the Taylor series around p = 0: 


i d k rr(j)/n \ 

H k (p,v)= V * + k = o,...,e-d, 


3=0 


where 


\Hk(p,<p)\ < cip £ d k+1 , p 6 [0,1], 

and ci > 0 does not depend on p 6 [0,1] and ip e <h. 

For £ — d > 0 and y > 0, we set 


(4.13) 


(4.14) 


Fk(y) := lr 


l-d-k 


Hk(p,<p) 


P 


l-d 


Hv) 


k =0 


E 


C,M)l := { 


£ — k — d^ (k — m)\ 

m =o v ’ 

£ ~ d H£~ d ~ m \0,y))r e - d - m 


E 

k m =0 


(l — d — m)\ 




l—d—k —1 

E 

3=0 


tfFW) 


//: 


j\(£-j-d-k) {£-d-k)\ 


- dp, k — 0, ...,£ — d, 

(4.15) 

J p£~d~k +1 

0 

(4.16) 

r k—m 

- if k — 0, — d — 1, 

(4.17) 

if k = £ — d, 


i-k)\ l nr )’ k = 0 ---’ e d ■ 

(4.18) 


Remark 4.2. 1. The integral in (I4.16P converges due to (14.14[) . 


2. The first formula in (14.171) is used only if £ — d > 1. 

3. In (I4.18P and below, we use the convention that the sum vanishes if the upper summa¬ 
tion limit is less than the lower limit. 

The following assertion follows from (1 1.1 1 1) and the definition of F k (y), F(y), and G k {y)- 
Lemma 4.3. Let £ — d > 0. 

1. The functions F k (y), F(y), and G k (y) are bounded near the origin. For each fixed 
they are continuous with respect to r > 0. 
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2. We have 


ft" E (ff) = F(x), x ez“\{0}. 

k=0 

Now, based on Lemma 14.21 in the case £ — d > 0, we find another preliminary asymptotic 
formula for v(x, t ) in terms of the above functions F k (y), F k (y), F(y), G k (y), and fl k (y). Now 
the asymptotic formula is extended to x — 0, but still contains an unknown grid function 
fl(rc), which will be determined in the next section. 

Lemma 4.4. Let £ — d > 0. Then, for any to > 0, integer K > max(M, £ — d), and t > to, 
the following hold. 

1. If x G Z d \ {0}, then, along with asymptotics (14.81) . we have 


i—d—1 


v(x,t) =5^ t 1 [d+k),e [P k [£ji) + G k (x)j + G e _ d (x) lnt 


k =0 


+ F*-d {p/i) + &( x ) + Hd+k)/t-i Fk (ji7l) + r A K ',x,t), 


(4.19) 


k=£-d -\-1 


t(d+k )/£-1 K Kt 1 / 1 J 


where 


2. If x = 0, then we have 


i-d 


&(y) = n(y) + ^n k (y). 


(4.20) 


k =0 


i-d -1 


v (0, t) — ^ + H e-d{ 0) Int 


k =o 


K 


+uj+ 


k=£-d -\-1 


1 £H k { 0) 
t (d+k)/i -1 i-d-k 


(4.21) 


+ r v (A';0,t), 


where u£ C does not depend on t 0 and t and F k ( 0) are defined in (14.7(1 . 


In both cases, r v (K-,x,t) satisfies estimate (14.9p . 

Proof. First, we fix x G Z d \{0}. Consider the terms in (14. 8 p corresponding to k — 0,..., £—d. 
As we will see below, they (in general) tend to infinity as f —>■ oo. For each m G {0,..., £ — d}, 
we have 


f l-(m+d)/£p ( X 4 

m \tw) 


= £r 


l—d—m 


H m (p, ip) 




p 


—d—m+l 


dp + £r 


l—d—m 




P 


fi-d—m -\-1 


dp. 


(4.22) 
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Consider the first integral in (14.221) . Using (I4.13j) . we represent it as follows: 




rt 1 /^ 

which yields 

i 


P 


1—d—m+l 


d P= 


£ d m 


P j 


-\-d-\-m—£—l 


dp + 


Hm{P,V) 


3 =0 


rt 


rt 


P J 


d—m+1 


dp , 


rt 


H (n in\ ^ C [ m 1 Vn rr(f-d-m) 

n m{P, <P) _ \^ Jim fu, _ ,l-(j+d+m)/^ , _ 

j\(£ — j — d — m) (.t — d — m)\ 


P 


l £—d—m +1 


R)™' i; (0,<p) In t 


3=0 


£—d—m—l 


- E 


h£\ o,ip) 


Tj(£—d—m) 
dl m 


(o,<p) 


j=0 


j\(£ — j ~ d — m) {£ — d — m)\ 


In r + 


Hm(Pi P) 


rt 


P‘ 


JL— d— m+1 


dp. 


Combining this with (14.22)1 . we obtain 

£ ^ 1 TJ'i.j') f \ 7 

/J /77 (0^) /'■ ^l-fd+m+j)/* 


(_£_) = t v 


3=0 

jj(£—d—m) i 

dim 


j\{£- j -d-m) 


(4.23) 


hit + n m (x) + ti-^/^m 


{£ — d — m )! 


KPPJ 


where F m (y ) are given by (14.151) and £l m (y) by (14.181) . 

Equalities (14.8ft and (14.231) imply (14.19ft with Gk(y) given by (14 .1 7ft and Cl(y) by (14.20ft . 
For x = 0, formula (14.21ft directly follows from integrating (14.101) with x — 0. □ 

The advantage of the representation (14.19ft is that it allows one to distinguish time- 
independent coefficients at positive powers of t in the asymptotic formula for v(x, t). Namely, 
the following assertion holds. 

Corollary 4.1. Let the assumptions of Lemma l4~4l hold. Then, for each x G Z d \ {0}, 

£—d—l 

v(x, t) — 7 1 -(d+A;)/£^ fc ( a; ) _j_ Gi-d(x) lnt + F(x) + Cl(x) + o(l) as t —y oo, (4.24) 

k=0 


where G^fx), F(x), and Cl(x) are given by (14.171) . (I4.16p . and (14.201) . respectively, and o(l) 
depends on x. 

Proof. The assertion follows from (j4.19j) and Lemma 14.31 (part [2]) • □ 

Remark 4.3. Asymptotic formulas (14.8ft and (14.19ft are uniform with respect to x G Z d \{0}, 
in the sense that the constant R V (K , t 0 ) in the estimate (14.9ft does not depend on x G Z d \{0}. 
Unlike those formulas, the function o(l) tends to 0 as t —» oo in relation (|4.24[) uniformly 
with respect to x 0 from any compact set, but not uniformly in Z d \ {0}. 
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However, in the next subsection, we will find an asymptotics analogous to (I4.24p directly 
from the integral representation (I2.23p of the Green function. In this new asymptotics, the 
time-independent term will be explicitly present and thus will coincide with F(x) + f2(x) 
from (I4.24p . Hence, we will find 0(a?) and thus 12(x) and “finalize” the (uniform) asymptotic 
formulas f!4.8p and (14.191) . 


5 Asymptotics of the second Green function v £ (x, t) 

5.1 Case t — d < — 1 

First, we complete Lemma [4.21 bv finding f2(x) in the case £ — d < —1. 

Theorem 5.1. Let i — d< —1. Then, for any e > 0, to > 0, integer K > M, and all x G 
and t > t 0 e e , 


K 

E 

k=M 


T F t (tTJi) + T UK;x,t), 


(5.1) 


f(k+d)/e -1 K \f_i!i) 


where 


Fk(y) are given by (14.5j) . 


Q(y) : = 


JyO 


d6, 


|r £ v (K',x,t)\ < 


(2v r) d JA{6) 
£ K+1 R v (K, t 0 ) 


(5.2) 


t(K+d+i)/e -i ’ 

and R v (K,t 0 ) > 0 does not depend on £ > 0, x e and t > t 0 £ e . 

Proof. Due to (I2.24p . it suffices to prove the theorem for e = 1 and v(x, t) = v^a^t). 
Since £ — d < 0, it follows from (14.7p that 


t (k+d)/e-i Fk ( t i/r) 
Therefore, by Lemma [4.21 


—» 0 as t —y oo, k = 0,1,2, 


f2(x) = lim v(x,t), x G if. 


t—> OO 


On the other hand, representation (I2.23P with £ — 1, estimate (I2.19p . and Lebesgue’s 
dominated convergence theorem imply that 


Jyd 

lim v(x, t) = . , , . . 

t^°° (27r) d J Rn A(9 ) 


dO , x G if, 


where the integral is absolutely convergent for £ — d < —1 due to (12.191) . The last two 
equalities, combined with Lemma [4.21 and Remark 14.11 complete the proof. □ 
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As in Corollaries 13. 1 1 and 13.21 Theorem 15. II and Lemma 14.11 item[3l imply that the higher 
approximation order M (see Definition 12. 3 p we choose, the faster the second Green function 
converges to the main term in its asymptotic representation, as t —» oo or as £ —> 0. 


5.2 Case i - d > 0 


First, we take £ = 1 and fix an arbitrary x G Z d . Using representation (I2.23P with e = 1, we 
write 

(5.3) 


v(x,t) = + v 2 (x,t) -v 3 (x,t)), 


VlOM) : = 


v 3 (x,f) : = 


l — g ~tA{6) 


e lx6 d9, v 2 (x,t) 


„ix9 


/ B 


-i /t 


m 

p—tA{6) 


’ Rtt\B 1/e A (0) 


de, 


e ixe de , 


iRn\B t _ 1/e A (9) 

where B\ := {9 G R d : \9\ < A} and R n is the cube defined by (11.lip . 

1. Estimate of vi(x, t). Making the change of variables 6 = we obtain 


Vl 


(t) = t~ d ' 1 


f 1 _ e -tAft- 1 /*® 

l Bl ~ 4|FnT 






(5.4) 


(5.5) 


Using the Taylor expansion for e lx6 and formula (I2.15p . we represent 


ixO ^ ^ 

x ) +n. 


(5.6) 


where Q-i(9,x) = ———, Qk-e{9,x ) are homogeneous functions of degree k — i with respect 
M(fc') 

to 6 representable as ratios of homogeneous polynomials, and there is C\ = Ci(x) > 0 such 
that 

\Q(0,x)\< Cl \6\- d+ \ 9 G R n . (5.7) 

Relations (15.6p and (15.7p yield 

= ixt e ~ d 


= X/ {k ti/e< ^ k ~ e ^ x ^ + ^ t 


where 


I Q(t l, ^ 1 x)\ < dt {d 1)/£ |£| d+1 , £ G R nt in, t > t 0 . 

On the other hand, due to (13.8p and (13.111) . 

i - = 1 - e-4«> - e"4(« ( r m/, V m+ ,K) + P(£,<) ) . 


(5.8) 


(5.9) 


(5.10) 


\m=M 
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where 


\P(S,t)\ < c 2 t-( K+1)/t \£\ K+1+e , £eB u t> to- (5.11) 

Combining (15.5|) with fl5.8pMI5.llD and recalling that K > £ — d, we obtain 


k =0 


l-d-M K 


B i 


- t 1 ( fc + m + o! )/^ f e A ®Q k _ e (£,x)V m+ t{€) d£ + o(l) as t —>■ oo. 

_n _A /f 


k =0 m=M 


B i 


Therefore, explicitly writing a coefficient at the zeroth power of t, we have 
vi(z,f) = v*(z,t) + f (1 - e - " 4 ®) Q_ d (£,z) 


Bi 


l-d-M 


(5.12) 


- X e AI -^Q k _ e (£,x)V 2 e-d-k(l)d£ +o(l) as 1 -f oo, 


k=0 


B i 


where v*(x, t) is a linear combination of the positive powers y 1 —(fc+rf)/^ ^ — o,..., £ — d — 1) 
with ^-dependent coefficients. 

2. Estimate of v 2 (x,t). We write Qk-i[0,x) (defined in (15.61) 1 in the spherical coordi¬ 
nates (r, p) with respect to 6: 


Qk-lifi,x) =r k l Q k -t(<P,x) 


(5.13) 


and note that Q k -i{p,x) are infinitely differentiable with respect to p and x. Recall the 
function r 7T (p) deffiied in (1 1.12 [) . Then, using (15.41) . (15.61) . and (15 .1 31) . we obtain 


e-d 




Q(9, x) dO. 


v 2 (x,t) = y; J Q k _ £ (p,x)J(p)dp J r d+k e 1 dr + 

k=0 $ t" 1 /* Rn\B t _ 1/t 

Since Q(9,x ) is integrable at the origin due to (15.71) . we calculate (as t —> oo) 

Inf 


Vo(x,f) = V 2 (x, t) + 


£ 


e-d -i 


^ t-d-k 


k =0 


Q- d (p,x)J(p) dp 

{r w (p)) d+k ~ e Q k - £ (p,x)J(p)dp + f Q(9, x) d9 (5.14) 


Rn 


In (r 7r (p))Q- d (p,x)J(p) dp + o(l), 
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where v^x, t) is a linear combination of the positive powers ft-fa+d)/# (k — 0,..., i — d — 1) 
with x-dependent coefficients. 

3. Estimate of v 3 (x, t). Making the change of variables 9 = we obtain from (j5.4j) 


v 3 (t) = t 


- +~d/t 




e ixt d£ + t~ d/e 




y ixt 


A1/£ & 


'\ t l/L\ B l A ( f ^0 

where L is defined by (13.71b Therefore, using (I2.19p . we obtain similarly to (I3.22p 




v 3 (i) = t 


- +-d/l 


R t i/L\Bi A (* 1/e £) 


z lxt 1/e ^ d £ + o(l) as £ —>■ oo. 


(5.15) 


Using (I5.15p . (13.8p . (13. lip , and (15.81) . we have 


A' 


v 3 (x,t) = t 


- 




R ntl/L\ Bl 


i+ £ r”*/'p m+ ,® + PK,i) 


m=M 


't-d 


E* (fc ^Qk-tfax) + Q(t v ^,x) df + o(l). 


vfc=0 


Taking into account estimates (13.121) and (15.9p and recalling that K > t — d, we obtain 


K 


t-d 


v 3 (x,t) = 

= v* 3 (x,t) + 


e ~ A ( 1+ E t ~ m/t 'P»+t(0 ^-^Q k . e ^ x ) d£ + o(l) 

m=M ) \k =0 / 

t-d-M „ 

Q- d (t„x)d£+ E / e~ A ^Q k _e^,x)V 2 t-d-k(0 d ^ + °( 1 ) ; 

;—n J 


-M), 


k =0 


(5.16) 

where v 3 (x, t) is a linear combination of the positive powers (A; = 0, — d — 1) 

with x-dependent coefficients. 

4. Estimate of v(x, t). Combining (15.3p with (15.121) . (15. 14p . (15.161) . we have for each 
x G 7L d 

In t I* 

J Q-d(v,x)J(v)d<p + S{x) + o{l), (5.P7) 


v(x, t) = v*(x,f) + 


where v*(x, t) is a linear combination of the positive powers t 1 ( fe + d m [k — 0, d — 1) 
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with ^-dependent coefficients and 


S{») ■■= 


(2ir) 


(1 - e-^0) Q_ d ((, y)d£- e~ A «>Q_ d ((, y) % 


y B i 


l-d-M 


Q(9,y)d9- y le A{0 Q k - e (£,y)V 2 <>-d-k(Q d£ 


Rtt 

e-d-i 


k =0 


S e-d-k I ( r ^)) d+k e Qk-£(v,y)J(v)d<p 


k =0 


ln ( r *(¥’))<5-i(»’> »)■%>) i* 3 


(5.18) 


Now we are in a position to formulate our main result in the case £ — d > 0. 

Theorem 5.2. Let £ — d > 0. Then, for any e > 0, to > 0, integer K > max(M, i — d), and 
all t > to£ l , the following hold. 

1. //i6Rj\ {0} ; then 


x 

t T / £ 


i-d 


v-w,i) = t'-«‘F 0 -U + v AI-M'F, + £ '-“n - 


if 

E 


£ 


k=M 

k 


X 

t~T 


t{k+d)/i- 


X 


~ Fk {jrji) 


k=ma,y;(M,£—d-\-l) 

where F k (y ) are gwen by (14. 4 p and (14.5H . respectively, 

l-d 


(Ky)=S(y)-F(y)~Yi 


k =0 

S^y) zs given by (15.18f) . -F(y) by (14. 161) . and h2 fc (g) by (I4.18|) . 
2. If x = 0, then 

i-d -1 

~kj.l—(d+k)/i 

l-d-k 

e k £H k ( 0 ) 


(5.19) 


£—d—\ 

V £ (0 ,t)= y £ k t i-(d+Wt eHk(0) + e <-d(ff g _ d(0 ) Illt _ ^_ d (0) In £ + ff(0)) 


fc=o 
if 

E 


t(d+k)/i-l £ _ d — k 


+ 4(X;M), 


k=£-d -\-1 

where S'(O) is given by (j5. 18[) with y — 0. 


(5.20) 
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In both cases, 


K(K]x,t) | < 


£ K+1 R v (K, t 0 ) 


f(K+d+i)/e-i 

and R v (K,t 0 ) > 0 does not depend on e > 0, x G Mf, and t > t^e 1 . 

Proof. The result follows from Lemmas 14.21 and 14.41 Corollary 14.11 and formula (15.1 7ft . 


□ 


Remark 5.1. Using Theorems 13.11 and 15.21 one can easily obtain asymptotic formulas for 
the first-order “spatial derivatives” £ _1 5 £ ,fc±u £ (;r, t) and £~ l 8 Ei k±v e (x,t), k = 1 and 

for their higher-order analogues. To do so, one should use Taylor’s expansions around x/t 1 ^ 
of the functions entering formulas (13.111 and (I5.19p . We will not provide this result in full 
generality, but rather present it for a one-dimensional example in Section 16.51 


6 Higher order approximations of the ID Laplacian 

In this section, we illustrate our general results by applying them to the one-dimensional 
discrete Laplace operator. Throughout this section, we fix N e N. In what follows, we will 
not explicitly indicate the dependence of constants, functions, etc. on N. 


6.1 List of constants and functions 


For asymptotic formulas below, we will explicitly find all the coefficients. In this subsection, 
we collect the constants and the functions which are needed for that. 

For integer J > 0, set 


bn ■ — 


2 (— 1 ) 


n+l 


N 


On m. • 


Vfl Z/2U+ 1 ) 

(2(n+l))!^ * 

E b h ■ ■ ■ ■ ■ 


ll~\ - \~lm—n 


^2 (n+m) 


ml 


P»K):= E ■ 

l<m<n/N 

dju ■= E (- 1 )' 

l<m<min(n/N,J) 


J 


m 


n > N, 

(6.1) 

n > N, m > 1, 

(6.2) 

n > N, 

(6.3) 

n > N, 

(6.4) 


Qjn(Z):= 0, n — N,..., 2N — 1, (6.5) 

Qjn(Z)--= E n > 2N , (6.6) 

h,h>N, 

h+l2=n 

1Z Jn (0 := (-1)' J e J (d Jn e n + V n (£) + Qjn(O), n > N. (6.7) 

As before, we use the convention that the sum vanishes if the set over which the summation 
is taken is empty. 
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6.2 Setting 

Let u £ (a;), s 6 Ij, be a grid function. We consider the difference operator 


N 


A £ = — A e , A e u(a;) := e 2 a p (u(x — ev) — 2u(x) + u(x + eu)) , (6.8) 


V=1 


with the coefficients a u satisfying the equalities 


N N 


a y v 2 = 1, a u v 2m = 0, m = 2,..., N. 


(6.9) 


The matrix of system (16.9ft is the Vandermonde matrix, which guarantees the existence and 
uniqueness of a v . Thus, throughout this section, £ = 2 and d — 1. 

The following lemma shows that Condition 12.21 holds with M = 2 N. 

Lemma 6.1. The operator A e is the (2 N)th order approximation of the second derivative 
(the one-dimensional Laplacian ) in the sense of Definition al 73l 

Proof. For any smooth function u(y) and e > 0, equalities (j6.9[) imply 

N 

s~ 2 ^2 a v{ u (y - £u ) - 2 u (y ) + u (y + £i> )) 

v=l 



(■ u [2N+2) (y - £ lv ev) + u (2N+2) (y + 6^)) 


where 



v=l 


□ 


Due to Definition 12.41 the symbol of the operator A e = — A e is given by 


N 



( 6 . 10 ) 


The following lemma shows that Condition 12.31 holds. The proof is given in Appendix lAl 


Lemma 6.2. Let A(6) be given by (16.101) . Then 
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1. W(0) = A'(it) = 0 and A'{6) > 0 for all 6 G (0,7r); 

2. A(O) = 0 and A(9) > 0 for all 6 G (0, 7r] ; 

K 

3. A(9) = 9 2 — b n 9 2n+2 + 0(9 2 ( R+2 l), where b n are given by (16.ip . 

n=N 

From now on, we confine ourselves by the case e — 1. The case e > 0 can be easily 
obtained by applying the scaling rule (|2.24p . Let S(x) := S 1 (x)^ x G Z, be the grid delta- 
function dehned in (1 1 . 2 j) with e — 1 . 

Definition 6.1. Let e — 1. We call the function u(x, t), x G Z, t > 0, the first discrete 
Green function if u(-, t) is a rapidly decreasing grid function for all t > 0, u(a;, •) G C^O, oo) 
for all iGZ, and 

f u(x,t) — Axufx, t) = 0, x G Z, t > 0, 

| u(a;, 0) = <5(x), iGZ. 

We call the function v(x,t), igZ, t > 0, the second discrete Green function if v(-,i) is 
a rapidly decreasing grid function for all t > 0, v(;r, •) G C^O, oo) for all iGZ, and 

J v(x,t) — A 1 v(x,t) = 3(x), iGZ, t > 0, 

| v(x, 0) = 0, iGZ. 


6.3 Asymptotics of the first Green function u(x,t) 

We begin with asymptotic formulas for the first Green function u(x, t). Consider the function 


h(y) ■■= 



( 6 . 11 ) 


which coincides with H(y) dehned in (12.71) . Note that —=h(y) is the Green function of the 

Vt 

operator u t (y,t ) -u yy (y,t) (cf. (EH). 

Theorem 6.1. For any to > 0, integer J > 0 and Ii± > N, and all x G Z and t > to, we 
have 


d J u(x, t) 
dt J 


GVt 


h 


JO 


x 

71 


Ri 1 

y 1 

“ t n+J Vt 

n=N 


h 


Jn 


X 

71 


+ r u (J,Ki]x,t), 


( 6 . 12 ) 


where 


hj 0 (x) := 


d 2J h( 


x 


dx 2J ’ 
r u (J,Ki,x,t)\ < 


hjni,•£) ■ F-Jn 

Ru(J, Ki,t 0 


4 ) ,lW - 


(6.13) 


t J+Kl+i \/t ’ 

n jn(0 are polynomials given by (16. 7$) , and i? u ( J, K 1 , t 0 ) > 0 does not depend on t > t 0 and 
iGZ. 
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Proof. Due to Lemma 16.21 part [21 the Taylor expansion of the symbol A{6) contains only 
even powers of 6. This fact and Theorem 13.11 with M = 2N, K = 2 Ki, and k = 2 n imply 
the result. □ 


6.4 Asymptotics of the second Green function v(x,t) 


For y > 0, we introduce the functions 


fo(y) 



p 2 h(p)dp, 


(6.14) 


fn(y) = -Jn [ P 2n 2h 0 n(p) dp, 71=1,2,.... (6.15) 

y Jo 

Note that 

fo(y) = h(y), y> 0, (6.16) 

where h(v) is riven bv (16. lip. Furthermore. fn(v) is positive, real analvtic foiPl y > 0 (e.g., 
use (16.161) ). and vanishes at infinity, together with all its derivatives. 

Using (16.lip and (16. 13p . we see that f n (y ) are real analytic for y > 0 and vanish at 
inhnity, together with all their derivatives. 

Theorem 6.2. For any t 0 > 0, integer K\ > N, and all x 6 Z and t > to, we have 


v{x, t) = Vtfo + D(x) + fn + x ^ °> 


v(x,t) = v(—x, t), 
where 


. 1 Iffcos(xf) 1 \ 2 \ x 

= tA.I hbrU + 2’ 


(6.17) 


X < —1, 


(6.18) 


fo(y), fi(y),---, Ik, (y) are given by (j 6.14j) , ( |6.15[) , 


r v (Ki,x,t)\ < 


R V {K i, t 0 ) 

t^Tt 


and R v (Ki, to) > 0 does not depend on t >t 0 and x > 0. 
Moreover, if N = 1, then f2(x) = 0 for all integer x > 0. 


Proof. Since v(x,t) = v(—x,t), x G Z, it remains to prove the first formula in (16. 171) . Due 
to Lemma [6.21 part [21 the Taylor expansion of the symbol A(9) contains only even powers 

1 Analyticity at y = 0 is understood in the sense that f{y) can be represented as the Taylor series with 
respect to the powers of y which converges to f(y) in a right-hand neighborhood of y = 0. 
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of 6. Therefore, H k (y ) = H ok (y ) = 0 for odd k. Hence, F k (y) = 0 for odd k , where F k (y ) 
are given by (14. 5p . On the other hand, relations Cj6.14j) . (16.151) and (14.4|) . (14. 5 1) . imply that 


fo(y) = F o{y), fn(y ) = F 2 n (y), n>i. 


(6.19) 


Thus, applying Theorem 15.21 with e — 1, M — 2N, K = 2/0, and k = 2 n, we obtain the 
hrst formula in (16.171) with f2(;r) given by 


fl(x) = — 
V ; 2t r 

Note that 




P 2 


a^dp-m 


( 6 . 20 ) 


kp) - m 


r 


dp — 2 


Kp) 


p z 


dp — h( 0) I = 


7T 


e _p2 / 4 - 1 1 

-5- d P = o> 

p - 2 


where the last equality can be proved, e.g., by introducing a parameter a as follows: e~ ap2 F 
and differentiating with respect to it. Hence, (I6.20P yields (16.18p . 

Finally, if AT = 1, then fl(a;) = 0 for all x G 7L by Lemma lB.ll □ 


6.5 Asymptotics of the gradient of the Green functions 


Next, we derive asymptotic formulas for Vv(i, t ) := v(x + 1, t) — v(x, t) and its time deriva¬ 
tives, see Remark 15.11 We introduce the function 


y{x) ■= foK), x > 0 , 


( 6 . 21 ) 


where f 0 (x ) is given by (I6.14p . Note that g(x) is negative, real analytic for x > 0, and 
vanishes at infinity, together with all its derivatives. 

Furthermore, for n > 2, we set 

9n(x) := ~Zvj fo n+1 \ x ) + x>0, (6.22) 

2s+j — l=n 


where f s (x) are given by (I6.15p . Note that f s (x) are real analytic for x > 0 and vanish at 
infinity, together with all their derivatives. Therefore, expanding the functions f s around 
x/\/i in (I6.17p . we arrive at the following theorem. 


Theorem 6.3. For any t 0 > 0, integer K\ > 2, and all x E Z and t > t 0 , we have 


Vv(x, t) 


VH(x) + g 





Vv(x, t) 


+ fl ^J2 9n + U(K 1]X , t ), 

= -Vv(-(x + l),t), 


x > 0 , 

x < — 1 . 


(6.23) 
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where £l(x) is given by (16.181) . g 2 (x),... ,g Kl (x) are given by (16.221) . 


r v (Ki',x,t)\ < 


Rv(Ki, to) 

t( K i + l)/2 ’ 


and, R v (Ki,to ) > 0 does not depend ont>t 0 and x > 0. 
Moreover, if N = 1, then VO(a;) = 0 for all integer x > 0. 

Finally, for n > 2, we set 


gjn(x) 


(n + 1 )! 


h 


(n+l) 

JO 


(x) + 


£ ^ 

s>N,j> 1 , J 
2s+j —l=n 


(i), 

Js 1 




x 


> 0 , 


(6.24) 


where h J0 (y ) and h Js (y) are given by (16.131) . Note that h J0 (y ) and h Js (y) are real analytic 
for x > 0 and vanish at infinity, together with all their derivatives. Therefore, expanding 
the functions hjo and hj s around x/\/t in (16.121) . we arrive at the following theorem. 


Theorem 6.4. For any t 0 > 0, integer J > 0 and Ki > 2, and all x E 7L and t > t 0 , we 
have 

d J ^(x,t) _ J_ {2J+1) 1 h(2J+2) (x) 

dM ~ t^ h [X)+ 2t^Vi [ ) 

K x 


£ 

n =2 


1 I X 

fJTTF^ 9jn 


d J X7 u(x, t) d J Vu(— (x + 1), t) 


dt J dt J 

where gj 2 (x),... ,gjK x (x) are given by (16.241) : furthermore, 


+ r u (J,Ki,x,t), x>0, 
x < 0. 


(6.25) 


r u (J,Ki,x,t)\ < 


Ru(J, Ki, to) 

t J+l+(I< 1 +l)/2 ’ 


and R u ( J, A\, t 0 ) > 0 does not depend on t >t 0 and x > 0. 


A Proof of Lemma 16.21 (Ellipticity) 


Let us prove assertion 1. 

Step 1. Let T u (y) and U v (y) be the Chebyshev polynomials of the first and the second 
kind, respectively. Using the relations 


cos u9 


7)) (cos#), 


sin u9 
sin# 


LU_i(cos #), 


vll v -x{y) = %(y), u — 1,2,..., 


we obtain 


N 

A'(9) = 2 sin # ajTf (cos #). 

v=\ 
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Thus, it suffices to prove that 


N 


®{y) '■= ^avViy) > 0 , ye [-1,1). 


u= 1 


Step 2. Suppose we have proved that 

(_l)i- 1 jB (j- 1 )(i)> 0> j — 1,N. 

Since T v (y) is a polynomial of degree v, it follows from (1A.2j) that 


(- 1 ) 


N 


_ l d N ~ l B(y) , 1) 


dy 


N 


^ = (- 1 ) 


dy 


JV 


^> 0 , ye [ 0 , 1 ]. 


(A.l) 


(A. 2 ) 


(A.3) 


Therefore, (—1) 


d N ~ 2 B(y) 

N ~ 2 7 is decreasing and, due to (1A.2|) . is positive for y e [—1,1]. 


dy 


N—2 


Hence, (—1) 


N _ 3 d N -*B(y ) : 


dy 


N—3 


is decreasing and, due to (1A.2I) . is also positive for y e [—1,1]. 


Continuing by induction, we conclude that B(y) is positive for y G [—1,1]. 
Step 3. It remains to prove (IA.2D . i.e., 


N 

(-1 r 1 5^ a v r} j \ 1) >0, j — 1,... ,N. (A.4) 

v=\ 

Set 

u{6) := 008 ( 6 *), v(y) : = arccos (y). 

Then T v (y) = u(vv(y)). In what follows, we will use the following representations for a = 

0 , 1 , 2 ,...: 

OO 

«<“>(») = (-i)“(i - y )- a+1/2 - </)", (A.5) 

n=0 

where A an > 0. This series and all the series below converge in a neighborhood of y — 1. To 
prove (1A.5I) . one can write arccos(y) = f*( 1 — ^) _ 1 A (1 + z)~ lB dz and expand (1 + z)~ x ^ 2 
into the Taylor series around z = 1. 

We fix some j G {1,..., IV}. From now on, we will not explicitly indicate the dependence 
of emerging coefficients on j. Since T u (y) = u{yv{y)), we see that 


% b) (y) = ^ (fe) Mz/)) B h-ik v{h) (y) ■ ■■■■ v{lk \y), (a. 6 ) 

k=1 1 <ll<-<lk, 

^lH- \-lk=j 

where Bi lmmm i k > 0. Therefore, using (1A.5D . we have 

j oo 

Ti>\y) = (-U V V C u W>(Mv))( 1 - y)‘- i+k/2 , (A. 7 ) 

k=1 1=0 
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where C k i > 0. 

Expanding u^(uv(y)) and the powers of v(y) (using (IA.5[) with a = 0), we obtain 


u {t \vv(y)) = V D km v m v m (y), 

m =0 

oo 

v m (y) = (l -y) m/2 Y, E ^ 1 -y) n ’ 


(A.8) 


n =0 


where D km e M, D km — 0 for k + m odd, and E mn > 0. Therefore, 


2 j—k oo 


U 


(fc) 


Ms/)) = EE Dkm Emn v m (l - y) n+m/ 2 + (1 - y) j - k,2 UUy ), (A.9) 


m=0 n=0 

where U k „(y) is continuous at y — 1 and 

f4i/(l) = 0. 

Now we substitute HA.9j> into (lA.7j) and take into account f)A. 10j> : 

j 2 j—k oo 

t?\v) = (-it E E E C k iD km E mn U k+m {l - y )l+n-3+{k+m)/2 + 

k =1 m=0 l,n=0 


where U v {y) is continuous at y — 1 and 


tE(i) = o. 


(A.10) 


(A.11) 


(A-12) 


Finally, we can calculate the sum in the left-hand side in f1A.4[) . To do so, we note the 
following: 

(a) in (lA.lip . there are only terms with is 2 , is 4 , ..., z/ 2j (since D km = 0 for k + m odd); 

(b) in (1A.11D . the terms with negative powers of (1 — y) cancel because Tv\y) is a poly¬ 
nomial (hence, contains only nonnegative powers of (1 — y))\ 

n 

(c) after summation with respect to is, the terms with is 4 ,, v 2j will vanish, while a v v 2 

V=\ 

will yield 1, due to flOlk 

(d) the terms with positive powers of (1 — y) and U v (y) will vanish at y — 1 due to (1A.12I) . 
Thus, only the terms corresponding to ( k,m ) = (1,1) and ( k,m ) = (2,0) remain: 

N ( 

(-l) i_1 I^ a ^ ) (l) = G-l) J '“ 1 (-l) i U > 11 Z C u E ln + D 20 J2 C * E ° 


J 0n 


v= 1 


n,l> 0, n+l=j —1 


n,l> 0, n-\-l=j —1 


(A-13) 
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Since u{6) = cos 9, it follows from the first equality in (I A. 8ft that D u and D 2 o are the 
leading order coefficients in the Taylor expansions of — sin 9 and — cos 9, respectively. Hence, 
Du = D 2 o = —1 and (IA.13D yields 

N 

Y, Cl lEln+ Y C 2 l E °n> 0, 

v=l n,l> 0, n-\-l=j — 1 n,l> 0, n+l=j —1 

which completes the proof o^| (1A.4D and assertion 1 in the lemma. 

Assertion 2 follows from assertion 1. 

The Taylor expansion in assertion 3 follows from the Taylor expansions of cos u9 and 
relations (j6.9j) . 


B One identity 


In this section, we prove a result, which we need for Theorems 16.21 and 16.31 in the case N — 1. 

Lemma B.l. Let 


Q(x) 


J _ I [ ( cos 00 
27 T 1 J V 2 ^ - cos 0 




x e z, x > o. 


Then h2(a;) = 0 for all integer x > 0. 
Proof. Step 1. First, we show that 


We have 


where 


h2(x) — 0(a: + 1) = 0 for all integer x > 0. 


n(;r) - + 1 ) = I(x) — -, 


T( s J_ f cos(a;0 -cos((a;+ 1)0 
An J 1 — cos £ 


(B.l) 


2 The left-hand side in (IA.4I) can be also found, using the following observation. In the proof, we have 

shown that, in (IA.6I) . only the terms with u' and u" are relevant. Furthermore, in the expansion of u'(ygfy)) 
only the term —vg{y) is relevant, and in the expansion of u"{vg{y)) only the term —1 is relevant. Therefore, 
we would obtain the same result if we replaced u(9) by the function —0 2 /2 and deleted all the terms with 
negative powers of (1 — y) after the respective expansions in the end. Therefore, 


(-D 


'7 1 ^ a 
1^=1 


-r«(i) = (-i 


V=1 


v arccos y 




v =i 2 dyi 


- (arccos 2 y) 


7 = 1 


because —^(arccos 2 !/) = 2 v(y)v'(y) already has no negative powers of (1 — y) due to (IA.5I) . Interestingly, 
ay 3 

the left-hand side in (IA.4I) does not depend on N. 
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If x = 0, then, obviously, I(x) = 1/2. 

Assume that x > 1. Using the formula 

cos(x + 1)£ = cosx£ cos£ — sin x£ sin £, 


we obtain 


I(x) 


1 

47r 


7r 



sin x£ sin £ 
1 — cos £ 



7T 

/ cos x£ d£ 


— 7T 


— 7T 


1 

47T 


7T 



sin x£ sin £ 
1 — cos £ 


d£. 


— 7T 


(B-2) 


For x = 1, we have I{x) = 1/2. To prove (IB. ID . it remains to show that the right-hand side 
of (IB. 2D does not depend on x > 1. Using the formula sin(x + l)£ = sin x£ cos£ + cosx£ sin£, 
we obtain 


smfx + 1)£ — sm x£) sm£ cosx£ sin 2 £ 

; -----= - sm x£ sm £ H---— 

1 — cos £ 1 — cos £ 

= — sin x£ sin £ + cos x£ + cos x£ cos £. 


For x — 1, the integral of the right-hand side vanishes due to direct calculation, while for 
x > 2, it vanishes due to the orthogonality of the systems {sinx£}a; e z and {cosx£} xe z, 
respectively, in L 2 (—7r,7r). 

Step 2. It remains to show that 12(0) =0. We have 


where 


12 ( 0 ) 


— lim 

7T cr—>-0 


^Ji(cr) + / 2 (cx) 





h(a) 


hi?) 


7r 


a 


2(1 — cos£) 


1 [ 

4 J sin 2 (£/2) 


/• d£ 1 1 

J 


1 

2 


cot(a/2) 


—h O(cr), 
a 


(B.3) 


(B.4) 

(B.5) 


Combining (IB.3|h (IB.4[) . and (IB. 51) yields 12(0) = 0. 
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